Large Orbital Magnetic Moment and Coulomb Correlation effects in FeBr 2 
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We have performed an all-electron fully relativistic density functional calculation to study the 
magnetic properties of FeBr2. We show for the first time that the correlation effect enhances the 
contribution from orbital degrees of freedom of d electrons to the total magnetic moment on Fe 2 + 
as opposed to common notion of nearly total quenching of the orbital moment on Fe 2+ site. The 
insulating nature of the system is correctly predicted when the Hubbard parameter U is included. 
Energy bands around the gap are very narrow in width and originate from the localized Fe-3d 
orbitals, which indicates that FeBr2 is a typical example of the Mott insulator. 

PACS numbers: 74.25.Ha, 71.27. +a, 71.20.-b 



FeX 2 (X= CI, Br, I) are well known model systems for 
the study of antiferromagnetismB They show unusual 
metamagnetic behavior wherein they undergo a phase 
transition from the antiferromagnetic to a saturated 
paramagnetic phase as a function of external magnetic 
field and temperature. FeBr2 has a layered structure, 
where the ferromagnetic Fe layers are widely separated 
by non-magnetic halogen layers producing a quite weak 
interlayer antiferromagnetic interactions (T/v= 14.2 K). 
The long range antiferromagnetic ordering can be over- 
come by applying ari|-external magnetic field (~ 30,-kOe) 
parallel to ther-axis-la As a result, many theoreticalu and 
experimentaoj works on FeBr 2 have been concentrated 
on probing the axial magnetic phase diagram, to predict 
and understand anomalies near the antiferromagnetic to 
the paramagnetic phase boundary. 

In addition, the electronic structures of transition 
metal dihalide system are of interest because they show 
various types of interesting behavior due to. the strong 
correlation effect in the transition metal ionfl For exam- 
ple, they behave as a Mott insulator or charge transfer 
type insulator depending upon the relative size of charge 
transfer energy A and the intra-orbital correlation energy 
U of the (i-electrons. 

Although large amount of work is done to understand 
the anomalies in the phase diagram and the correlation 
effects, surprisingly no electronic structure is known from 
the first principles calculation. In this report, we study 
the electronic and magnetic properties of one of these 
systems (FeBr2) by combining a first principles method 
and the Hubbard on-site correlation term. 

Orbital magnetic moment plays a crucial role in de- 
termining many important effects in magnetic materi- 
als such as magnetic crystalline anisotropy, norucollinear 
magnetism, magneto-optical Kerr effect etc.. 13 It is a 
well known fact that crystal field interaction quenches 
the orbital magnetic moment in systems containing 3d- 
transition metal. Accordingly, it is assumed that because 
of strong crystalline field on Fe 2+ ions from the surround- 
ing Br^" , the orbital moment will also get. quenched in 
FeBr2.la Ropka, Michalski, and RadwanskiEl showed that 
the orbital moment is not quenched fully by the crystal 




FIG. 1: Crystal structure of FeBr2. Black and white spheres 
represent Fe and Br atoms, respectively. 



field through their quasi-atomic calculation. Ropka et al. 
ascribed the origin of the large orbital moment to spin- 
orbit coupling (SOC). Our study also demonstrates that 
orbital degrees of freedom has substantial contribution 
(of ~ 20% to the total magnetic moment). However, 
correlation effect in Fe 3d orbitals increases the orbital 
contribution to the total magnetic moment significantly, 
as discussed below. This is explicitly shown by incor- 
porating the Hubbard parameter U in the self-consistent 
calculation which also predict an insulating nature. 

FeBr2 is isomorphic to the hexagonal Cdl 2 structure 
(space group p3ml, D^ d , No. 164) at atmospheric pres- 
sure and low temperature. Figure [l] shows the crystal 
structure. i-Thc lattice parameters are a = 3.77 A and 
c = 6.22 AH where the c/a ratio (~ 1.65) is very close to 
the ideal value of y/8/3 ~ 1.63. The chemical unit cell 
has one Fe and two Br atoms. The Br atoms are bonded 



2 



tightly to the Fe atoms on either side by covalent bond to 
form a Br-Fe-Br trio. The Br-Fe-Br trios are separated 
by an empty layer between them, and bonding between 
the trios is a weak van der Waals type. The Fe ions are 
located at la positions whereas Br atoms are located at 
2c? positions in a unit cell. Each Fe ion is at the center 
of octahedron formed by the surrounding Br ions with a 
small distortion. As a result, a small trigonal component 
of crystal field due to this distortion is present apart from 
the cubic component. FeBr 2 is also taken as an example 
of dynamic Jahn- Teller system wherein strong coupling 
between lattice and spin excitations arepHianifested in 
temperature dependent Raman spectrumJlj 

We have carried out density functional calculations0 
on FeBr 2 using the full-potential linear muffin-tin orbital 
(FPLMTO)lia method with foil relativistic effects on the 
single-particle electron states and the nonrelativistic local 
density approximation (LDA) for the many-body effects 
in the exchanga-correlation energy with the Vosko-Wilk- 
Nussair form.Ej The strong correlation effect in the. Fe 
3c? site is taken into account by LDA+U methodo In 
the past, LDA+U method is frequently used to predict 
the occurrence of large orbital moment and insulating na- 
ture andis found to be a robust method for the correlated 
systems Jij'LJ The chemical unit cell is doubled (c = 12.44 
A) along c-axis to simulate the antiferromagnetic system 
with 2 Fe atoms and 4 Br atoms. Basis functions, elec- 
tron densities and potential were calculated without any 
shape approximation. The trial wave function, poten- 
tial and charge densities are expanded in combination of 
spherical harmonics (\ m ax — 6) inside the atomic muffin- 
tin sphere region and plane wave Fourier series in the 
interstitial region (E cut =94 Ry). 2 K-basis set (where n 
is kinetic energy of muffin-tin orbitals in interstitial re- 
gion) is used. In the Fe sphere, 4s, 4p and 3c! orbitals 
are taken as valence states and Zp orbitals are treated as 
semi-core state, whereas in the Br sphere, 3d orbitals are 
taken as a semicore state and 4s, 4p and 4g! orbitals are 
taken as valence states. The semicore states are treated 
in a separate energy window. The k-space integration 
is done with 30 k-points in irreducible part of Brillouin 
Zone(IBZ)(144 k-points in full Brillouin Zone(FBZ)) us- 



TABLE I: Direct energy gap(Eg), Indirect energy gap(E g ), 
orbital magnetic moment(/i;), spin magnetic moment(/i s ), 
and total magnetic moment fj, = fj,i+p, s for FeBr2. The energy 
gap is in eV and magnetic moments are in Bohr magneton. 
U= eV in LDA calculation. 
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ing the tetrahedron method.0 We have also carried out 
calculations with larger set of k-points (43 points in IBZ 
and 320 points in FBZ), but the total energies are found 
to be converged within 1 meV and magnetic moments 
change within 0.01 \ib- To bring about the orbital contri- 
bution to total magnetic moment, spin-orbit interaction 
is included in the crystal Hamiltonian. The calculations 
are done at experimental lattice constant. 

Due to mean field character of LDA, it incorrectly de- 
scribes the ground state of many strongly correlated mar-, 
tcrials, for example, the transition metal oxide systems.tLSl 
In these materials, the d-orbitals are well localized and 
retain the strong atomic-like character. For a good de- 
scription of strong on-site correlation effect between elec- 
trons in the d-shell, the LDA+U method is widely used. 
This method identifies these orbitals as correlated states. 
The LDA+U variational total energy functional takes the 
form 



E tot [p, n] 



E LDA [p] 



E u \h] 



(1) 



where the first term is the usual local spin density func- 
tional of the local electron spin density p a (r)(a =T,i). 
h a is a local orbital (d or /) occupation matrix. Fi U 
is an electron-electron interaction energy which depends 
on Slater integrals F°, F 2 and F 4 . E dc is a double count- 
ing term. The Slater integrals are in turn related to the 
screened Coulomb and exchange parameters U and J re- 
spectively: 



and 



F 4 + F 2 
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U = F°. 



(2) 



(3) 



"Reference [9] 
'Reference [8] 



For 3d electron systems it is found that F 2 /F 4 = 0.625. 
It is to be noted that the orbital polarization is automat- 
ically included in LDA+U method, as the orbital depen- 
dent effect comes from F 2 and F 4 . Since no experimental 
value of U in FeBr 2 system is available in the literature, 
we decided to take U and J values from experiments on 
other systems. It is to be noted that the experimental 
value of U in FeO systemlij is 5.7eV and 6.0 eV. In similar 
systems like antiferromagnetic transition metal oxidesjlil 
J is close to 0.9 eV and does not change much from this 
value, so we have taken J = 0.9 eV. To check the validity 
of our U values, we performed the self-consistent calcu- 
lations at four different values of U, namely, 5.7 eV, 6.0 
eV, 6.4 eV and 6.8 eV. We consistently obtain the in- 
sulating ground state in agreement with the experiment, 
although the energy gap changes by a small fraction. The 
magnetic moments are almost constant with respect to 
the U values. So, the four different values of U do not 
affect our conclusions. We present our results for U = 6 
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FIG. 2: Energy band structure of FeBr2 at U = 6.0 eV along 
high symmetry directions in IBZ of the hexagonal unit cell. 
The top of the valence band is set to energy zero. 

eV taken from the experimental value in FeO system if 
not mentioned explicitly. 

Table § shows energy gap and magnetic moment 
for FeBr2 at the two representative values of U. LDA 
and GGA both predict FeBi'2 to be metallic, while 
LDA+U predicts it to be insulating in agreement with 
experiments It is, therefore, crucial to take into account 
the correlation effects to obtain the correct electronic 
structure. It means that there is a strong intra-orbital 
correlation for Fe 2+ ion in FeBr2 and that this strong 
correlation is responsible for the insulating nature. 

Figure || shows dispersion curves for electronic states 
along high symmetry directions in IBZ of the hexagonal 
unit cell. The top of valence band is set to energy zero. 
Br-4s bands are not shown in the figure which is located 
at around -16 eV. It is seen from Fig. |^ that there is a in- 
direct energy gap between H and T with a gap size of 1.53 
eV. A direct gap of 1.87 eV is located at the middle of 
the symmetry line connecting M and T. No experimental 
value of the gap is available m the literature. However, 
since FeBr2 is yellow in color the gap would be around 
2 eV which may correspond to a transition across the 
direct gap. Our calculated gap values of 1.87 eV are in 
reasonable agreement with the conjectured one(see Ta- 
ble ||). We can see from Fig. [2|that the bands have little 
dispersion in the z-direction (see the bands along sym- 
metry lines T-A, L-M, and K-H) which is a characteristic 
feature in systems with layered structure like graphite or 
some high T c superconductors. The equienergy surfaces 
near the energy gap, therefore, have cylinder-like shapes. 

Density of states (DOS) of FeBr 2 in Fig. ||(a) shows 
the total DOS for FeBr 2 . Spin- up Fe-d states are located 
deep in energy between -10 eV and -6.5 eV, whereas spin- 
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FIG. 3: Density of states for FeBr 2 . (a) Total DOS (b) Par- 
tial DOS of Fe-3d states (c) Partial DOS of Br-4p ( solid line) 
and Fe-4s ( dotted line). In (b) and (c), positive and negative 
DOS represent spin-up and spin-down DOS, respectively. 



down Fe-d states are located at between -0.5 eV and 3 
eV as shown in Fig. 0(b). Fe-d orbitals form localized 
states with narrow band width where the main peaks are 
located -9.6 eV, -8.9 eV, -6.8 eV for spin-up states and 
-0.2 eV, 1.6 eV, 2.2 eV for spin-down states. A group 
of bonding orbitals with wide band width is located be- 
tween -6 eV and -2 eV which is the hybridization of Br- 
p and Fe-s bands( Fig. ||(c)). This wide band gives a 
covalent bonding within the Br-Fe-Br trio. The energy 
gap between the hybridized bonding orbitals of Br-p/Fe- 
s and the conduction band is 3.6 eV. The d-bandwidth, 
W, of 0.5 eV is found to be much smaller than U so that 
U/W > 1. The states near the band gap at conduction 
band minimum and valance band maximum are mainly 
from spin-down Fe-d states. This means that FeBr2 is a 
typical Mott insulator (materials in which the kinetic en- 
ergy gain is smaller than the Coulomb repulsion energy 
U and as a result electrons can hardly hop to the Fe-d 
orbitals). 

Fe 2+ ion in FeBr2 has six electrons in d orbitals. The 
nearly octahedral arrangement of Br~ surrounding Fe 2+ 
ions introduces crystal field which splits the d-orbital into 
tag (dxy, d yz , and d xy ) and e g (d x 2_ y 2 and d 3z 2_ r2 ) or- 
bitals. This splitting hinders the free rotation of the elec- 
trons and reduces its orbital moment. Both t2 S and e g 
orbitals of spin-up Fe-d bands are fully occupied while 
only one of t2 g orbitals is occupied in spin-down bands 
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to contribute 4 fis as a spin contribution to the total 
magnetic moment. This simple estimation of the mag- 
netic moment from DOS is consistent with the calculated 
spin contribution of 3.86 fiB to the total magnetic mo- 
ment. The Fe 2+ ion has L=2, S=2 and J=4 in free space. 
If we assume that the orbital moment is quenched as in 
usual 3c? transition metal magnetic system, the spin con- 
tribution to the total magnetic moment is also consistent 
with the simple -atomic picture. However, the neutron 
diffraction studya shows that the magnetic moment on 
Fe site is nearly 4.4 [Lq. Wilkinson et alB argue that this 
moment is close to that expected for divalent metallic 
Fe 2+ ions if the orbital contribution is taken to be zero. 
Our calculation shows that total magnetic moment /.i is 
split into orbital fj,i and spin parts fj, s with the orbital 
part contributing a significant amount. From Table |] it 
is seen that the orbital contribution to the total mag- 
netic moment is ~ 20% the total magnetic moment. It 
is interesting to note that the nearly octahedral arrange- 
ment of Br ions surrounding Fe ion is not able to totally 
quench the orbital degrees of freedom. The orbital and 
spin degrees of freedom are therefore not totally decou- 
pled (obeying Hund-'s third rule). Quasi-atomic calcula- 
tions by Ropka et aH also give the orbital moment close to 
our values, and the large orbital moment were speculated 
to be due to SOC. However, our calculations show that 
SOC is not enough to explain the large orbital contribu- 
tion as can be seen in table |. LDA calculation with SOC 
gives an orbital magnetic moment of 0.14 hb and total 
magnetic moment of 3.64 fj,s which is smaller than the 
experimental value of 4.4 /j,b- When we include the on- 
site correlation effect at the Fe site as well as SOC, we get 
about five times larger orbital magnetic moment of 0.66 
fiB- Total magnetic moment of 4.52 hb is in good agree- 
ment with experimental value within 3% error. Thus we 
can see that the important contribution to the large or- 
bital magnetic moment of FcBr2 is mostly due to the 



correlation effect of the 3d orbitals. Therefore, the corre- 
lation effect plays a very important role in determining 
the magnetic properties as well as the insulating ground 
state. In Ropka et a/.'s quasi-atomic model, the correla- 
tion effects are included by the correlated basis set with 
no explicit term representing on-site correlation effect in 
the Hamiltonian. However, in our calculation we include 
explicitly the correlation term in the Hamiltonian. 

We suggest-that, if an x-ray magnetic circular dichro- 
ism (XMCD)E2l experiment, which is a local probe for 
spin and orbital magnetic contribution to total magnetic 
moment, is performed on FeBr2, then the orbital contri- 
bution can easily be probed. Such experiments relate, via 
atomic sum rules, the dichroic intensities [ie., difference 
between right- and left-circularly polarized photons ab- 
sorption cross-sections) measured at the absorption edges 
of the constituent elements to the ground-state expecta- 
tion value of effective one-electron operators such as (L z ) 
and (5,). 

To conclude, we have shown for the first time that the 
correlation effect is important in describing the magnetic 
and electronic properties of FeBr2 which plays an im- 
portant role in enhancing the orbital magnetic moment 
and predicting the correct insulating nature. This indi- 
cates that Fe 2+ ion in FeBr2 retains to some extent its 
atomic-like character. FeBr2 can be classified as a typ- 
ical example of Mott insulator due to the character of 
the energy gap and the width of the bands around the 
gap. We expect that the present calculation will help un- 
derstand the electronic and magnetic properties of other 
iron halide systems. 
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